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Abstract 

We present an analytical approach to the calculation of the linewidth and lineshift of an atom 
or molecule in the near field of a structured dielectric surface. For soft surface corrugations with 
amplitude A/50, we find variations of the linewidth in the ten percent region. More strikingly, the 
shift of the molecular resonance can reach several natural linewidths. We demonstrate that the lateral 
resolution is of the order of the molecule-surface distance. We give a semiquantitative explanation of 
the outcome of our calculations that is based on simple intuitive models. 
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61.16.Ch - Surface structure 

78.66.-w - Optical properties of specific surfaces and microparticles 



1 Introduction 

It is well accepted that the natural linewidth of the 
excited state of an atom, as well as the exact value 
of its energy levels, are greatly influenced by quantum 
fluctuations. When an atom is confined in a geometry 
the existence of the boundary conditions for the elec- 
tromagnetic field results in modifications of the atomic 
radiative properties. There have been many theoreti- 
cal works on the calculation and interpretation of these 
effects using quantum electrodynamics Q. Much of 
the physics involved can be addressed, however, by re- 
placing a two-level atom by a classical dipole moment 
and treating its radiation in the presence of boundaries 
0] . Such an approach is quite successful in treating the 
modification of spontaneous emission due to the new en- 
vironment. The energy level shifts of the atomic states 
in the near field can also be described very well using 
this model || . One finds the well-known Lennard- Jones 
potential which is proportional to 1/z 3 . In the far field, 
however, the Casimir-Polder shift, as well as the ex- 
act numerical value of the oscillatory resonant coupling 
of the excited state can be obtained only from a fully 
quantum electrodynamic treatment fl Bj . 



From the experimental side many groups have tried 
to study various aspects of these phenomena in different 
systems. The first experimental evidence for the mod- 
ification of spontaneous emission was demonstrated in 
1970 by Drexhage §. Here a very thin layer of fluo- 
rescing ions were separated from the underlying surface 
by a thin spacer layer, and the emission lifetime was 
recorded for different spacings. This technique has been 
used extensively ever since due to its simplicity and its 
very high vertical spatial resolution ||. Direct exper- 
imental verification of the energy level shifts of atoms 
in confined geometries was also demonstrated success- 
fully more recently by performing high resolution spec- 
troscopy @, H Following the discovery of Surface 
Enhanced Raman Spectroscopy in the early eighties 
the more complex case of a molecule in the vicinity 
of rough surfaces attracted much attention. Several re- 
searchers have studied the emission properties of an en- 
semble of dye molecules on rough surfaces and gratings 
10, OJ [l2|, [TsJ . Very recently there have been also some 
efforts on the spectroscopy of atoms placed on a thin 
organic layer above a rough surface |l4j]. 

In this paper we treat the modification of the ra- 
diative properties of an atom or a molecule placed 
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very close to a surface with lateral optical contrast. 
Our work is mainly motivated by the recent progress 
in the field of Scanning Near-field Optical Microscopy 
(SNOM) which has opened the door to optical mi- 
croscopy and spectroscopy with lateral resolution be- 
yond the diffraction limit. In the most common SNOM 
configuration one arrives at this high resolution by ex- 
amining the sample in the near field of a sub- wavelength 
metallic aperture. Indeed, single molecules on a surface 
have been detected with this method, and it has been 
verified that the molecular lifetime is modified by the 
presence of the aperture fl5|| . A more elegant approach 
to SNOM uses the fluorescence of a single molecule 
as a probe ||l6[ . Here one can record the molecular 
emission intensity or alternatively the molecular life- 
time as a measure for the interaction of the molecule 
with the sample surface. Girard and coworkers have 
shown numerical calculations for the modification of 
the molecule's lifetime as it is positioned above a sam- 
ple with nanometric topographic features |l7], [l8|]. In 
the present paper we propose an analytical approach 
for this problem based on a perturbative method from 
scattering theory. Our approach is valid in the domain 
of soft surface corrugations where a complex surface 
geometry can be Fourier decomposed in terms of sinu- 
soidal surface gratings whose corrugation amplitude is 
small compared to the molecule-substrate distance. In 
addition to the modification of spontaneous emission 
we also consider the modifications in the molecular en- 
ergy level shifts. The latter is particularly interesting in 
view of the recent achievements in high resolution spec- 
troscopy of single molecules In Konstanz we are 
currently pursuing experiments which aim at the mea- 
surement of the energy level shifts of a single molecule 
in the vicinity of a surface [^(J . As we show in this pa- 
per one can take advantage of the extremely high lateral 
resolution in this system to perform a novel form of op- 
tical microscopy. 

2 Presentation of the model 

We are interested in the radiative properties of an atom 
or molecule (called 'molecule' in the following) at a po- 
sition r in an inhomogencous environment. In this sec- 
tion we first outline the description of the environment 
and then discuss the model taken for the molecule. 

2.1 Environment 

We consider the molecule to be placed in the vicinity 
of a solid substrate at a distance ranging from a few 
nanometers to a few optical wavelengths (sec Fig.Q). 
These distances are large compared to the dipole's di- 



mensions, ensuring its purely electromagnetic interac- 
tions with the surface. Moreover, it is appropriate to de- 
scribe the solid by a local dielectric function e(x; w), al- 
lowing the description of the electromagnetic field phe- 
nomenologically by Maxwell equations. 

The substrate surface is given by the equation x% = 
s{xi,X2) where x = (x±, X2, X3) are cartesian coordi- 
nates. The surface corrugation is characterized by two 
length scales: the typical height h of the vertical corru- 
gation and its lateral scale of variation a. In the case 
of a surface relief grating h would be the grating ampli- 
tude and a its period. In the region below the surface, 
X3 < s{x\ 1 X2), the dielectric function is equal to the 
squared index of refraction n 2 that is assumed to be 
real. The dipole is located in vacuum at the position 
r = (x, y, z). We explicitly retain its lateral coordinates 
x, y since we are interested in the lateral resolution ob- 
tained in the variations of the linewidth and lineshift. 
In order to simplify the formulas, we shall use the no- 
tations X = (x\,X2) and R = (x, y) for the lateral 
coordinates. 

2.2 Linewidth and frequency shift 

To begin with, let us focus on a system with two energy- 
levels, a ground state E g and an excited state E e con- 
nected by an electric dipole transition. We assume that 
the transition dipole is oriented parallel to the x j -axis 
(j = 1,2,3, linear polarization) and write T>j for the 
dipole matrix element. When this atom interacts with 
the electromagnetic field its energy levels get shifted by 
amounts SE g ^ e (r), and the excited state acquires a finite 
lifetime 1/Tj(r) where Tj(r) is the spontaneous emis- 
sion rate. Both linewidth and level shifts may be calcu- 
lated in second-order perturbation theory. One obtains 
the spontaneous emission rate [plf 

r j (r) = -^- J dr(vac|^ (+) (r,i+r)£;^ ) (r,t)|vac)e i ^ T 

— OO 

(1) 

where E^^r, t) are the positive and negative frequency 
parts of the electric field operator at the atom's posi- 
tion r, and oj cg — {E c — E g )/h is the atomic transi- 
tion frequency. At this point one often proceeds to 
a mode expansion of the electric field operator, and 
the linewidth (|l|) connects to squared mode function 
amplitudes. We take here a different route, follow- 
ing the response theory developed by Agarwal p2| , 
and Wiley and Sipe More specifically, we in- 

voke the the fluctuation-dissipation-theorem to con- 
nect the linewidth Tj to the classical Green function 
Gij(x, r; lu). This Green function describes the electric 
field E d ip(x)e- iwt + c.c. (the 'dipole field') created by 
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Figure 1: Geometry of the problem. The surface X3 = s(xi, X2) separates the dielectric below (index n) from vacuum 
above. Three different length scales are involved: the dipole's distance z from the mean surface, the vertical surface 
corrugation h and the lateral corrugation scale a. The transition wavelength A is not shown: typically, z, h, a < X 
in near-field microscopy. 



an oscillating point dipole de luJt + c.c. located at r: 

£di P , 4 (x) =^2G ij (x,r;u)d j . (2) 
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The fluctuation-dissipation-theorem allows one to ex- 
press the linewidth (Q) in the following form |22|, 



— Im G j:j (r,r;u eg ). 



(3) 



In the vicinity of an interface the electric field radiated 
by the dipole differs from that in free space: it contains, 
in addition to the well-known dipole field , the field 
reflected from the surface. We write this field in terms 



of a Green function G^(x, r;w )- Upon insertion into 
Eq.(|^), we find the environment-induced modification 
STj(r) of the linewidth, that now depends on the atom's 
position relative to its inhomogeneous environment. We 
stress that in the present approach, the linewidth is lin- 
early related to the electric field radiated by the dipole, 
and it is not necessary to compute squared field mode 
amplitudes which is a more difficult task in a complex 
geometry. 

Let us now turn to the shift of the atomic resonance 
frequency dcj eg j(r) = (SE c (r) — 5E g (r))/H in the pres- 
ence of an interface. By a calculation similar to the 
one for the linewidth, Fermi's Golden Rule yields the 
following result (obtained from Eq.(2.9) of Ref.pl): 



Slo c 



V 2 , 



Re Gjj(r,r;w g) + 



duj Im Gjj (r, r; u) 

2lT (Jj 0£ r + UJ 



(4) 



The first term has a similar form as Eq. (^) but involves the real part of the Green function. 

For an atom with more than two states one has to take into account allowed dipole transitions to other energy 
levels. Let us focus on the situation depicted in FigJ| where the excited level E e is the first level above the ground 
state E g . The linewidth Tj is then still given by the two-level expression ([!]), and the lineshift contains an additional 
contribution from higher-lying states Ey, 



Su othcl 



E 

b,E h >E e 



dw 
2^ 



Im Gjj(r, r; 



UJ 



l(bKlg)P 



(5) 



If the excited state decays to more than one lower-lying 
level, the decay rate Tj is a sum over several contribu- 
tions, each one of the form (||). 

It is instructive to compare the linewidth (^) and 
the frequency shift (^) obtained from quantum theory 
to the corresponding quantities for a classical harmonic 



oscillator. This model of the Lorentz atom is widely 
used in the optics community [|[ [|, |||, ||(| , and we can 
make contact with the work done there. The Lorentz 
atom is a harmonic oscillator driven by the local elec- 
tric field, i.e., the external field plus the dipole field (0). 
Since this field is proportional to the dipole moment 
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Figure 2: Sketch of a multilevel atom with ground and first excited state. 



itself, it shifts the resonance frequency and leads to a 
finite damping rate. 

For a linearly polarized dipole along the Xj-axiB posi- 
tioned in an inhomogeneous environment, the shift 5u>j 
and the damping rate Tj are obtained from a straight- 
forward calculation [^J. Normalizing to the free-space 
linewidth T^, one obtains for weak radiation damping 
(in SI units) 



i lj(r) 

2 r 



37T£(; 
1,3 



G#(r,r;o;o) (6) 



where fco = ^o/c is the vacuum wavenumber. As in the 
quantum mechanical calculation, the imaginary part of 
the classical Green function determines the linewidth. 
The fluorescence rate Tj may thus be computed clas- 
sically ]2l]| , and although we focus in the following on 
the classical dipole model, our results for the linewidth 
remain valid for a generic atom. As for the lineshift 
given by Eq.(||), the classical calculation only yields the 
first term of the quantum-mechanical results (Q, ^) , and 
the nonresonant frequency integrals of the Green func- 
tion are not accounted for. This implies a limited va- 
lidity of the classical model for frequency shift calcula- 
tions. In this paper, however, we are mainly interested 
in studying the variations in the radiative properties of 
an atom as a function of its lateral position very close 
to a structured surface. It is quite common that the 
atom-surface distances are much smaller than the rele- 
vant atomic transition wavelengths. In this regime, the 
full quantum-mechanical lineshift (^, ^) approximately 
yields the electrostatic result 



(e\£ 



(g\d 2 M 



ReGii(r,r; W = 0), (7) 



Note that this shift is again determined by a classical 
Green function, in this case at zero frequency The 
prefactor is different, however, from the classical dipole 



and involves the difference in size of the electronic wave 
function in the ground and excited states (a difference 
that cancels for a two-level atom [^8| ) . As far as the lat- 
eral resolution is concerned, we may therefore treat the 
case of a classical dipole and use Eq.(||) for simplicity. 
The exact magnitudes of the line shifts and linewidths 
for realistic molecules can be then easily calculated con- 
sidering the above-mentioned discussion. 

3 Reflected field calculation 
3.1 Outline 

Our task is now to compute the Green function above 
the substrate, i.e., the reflected field created by an oscil- 
lating point dipole. The basic equations are the macro- 
scopic Maxwell equations, given the dielectric function 
e(x;w) and the external current j(x, t) — — iui d8(x — 
r)e _I "°* + c.c., supplemented by the boundary condi- 
tions for the field at the surface. We make the following 
ansatz for the electric field in vacuum above the solid 



E(x,<) = (E fs (x)+E r (x)) 



(8) 



where E fs (x) is the dipole field in free space, and E r (x) 
is its environment-induced modification (the reflected 
field). The latter is source- free above the surface. The 
total field (||) is matched at the surface of the solid 
to a 'transmitted field' (source-free inside the solid). 
This matching determines the reflected and transmit- 
ted fields in terms of the free space dipole field. 

In general the boundary conditions are complicated, 
and exact solutions are only known for simple geome- 
tries. In order to proceed analytically, we resort to an 
approximate solution for a 'slightly corrugated surface', 
i.e., a vertical corrugation small compared to the sep- 
aration of the molecule from the surface. The precise 
validity of this approximation is discussed in SecJ|. The 



4 



boundary conditions may then be linearized around the 
mean value of the surface function and solved to first 
order in the surface corrugation. We thus end up with 
three terms for the linewidth 



r 3 -(r) = r 00 + *r§(«) + irJ(R,«) 



(9) 



The first term is the vacuum linewidth. The second 
comes from the reflection at the flat surface (zeroth or- 
der field), and the third from the first-order scattering 
off the surface corrugation. This last term contains in- 
formation about the lateral surface structure since it 
depends on the lateral coordinate R = (x, y). A re- 
lation similar to (||) also holds for the frequency shift 
<Wj(r). 

3.2 Approximate solution 

In order to find the reflected field to first order, we 
use the so-called 'method of small perturbations' well- 
known in light scattering from rough surfaces [30], 
[3l[ . With this method, one usually computes the re- 
flected field for an incident plane wave. We thus expand 
the free space dipole field into Fourier components and 
work out the reflected field for each component. This 
approach is similar to that of Rahman and Maradudin 
32 1 and of George and co-workers fiT 
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3.2.1 Fourier expansions 

The free space dipole field is well known and may be 
found from the following formula for the Green func- 
tion (in SI units) H, 



d , r , 2 A expifcp x-r 
a — ^ \-Oijko — : (10) 



with the wavenumber fco = lo/c. Using the notation 
K = (fci,^) for the lateral wave vector components 
and recalling the notations x = (X, X3), r ~ (R, z) 
etc. for the lateral and vertical coordinates, the Weyl 
expansion reads pl|: 



expifco|x — r| 
47T£o|x — r| 



d 2 K 



2tt) 2 2e fc 3 
x exp i[K • (X — R) + fc 3 \x 3 — z 



(11) 



Here, the vertical wave vector component fc 3 is defined 

by 

k 3 = Jkl - K2 (12) 



and the square root is chosen such that Re k 3 , Im k 3 > 
0. It is important to note that the expansion ( |Tl| ) con- 
tains both 'far field' and 'near field' contributions, cor- 
responding to lateral wave vectors with magnitude K 



smaller and larger than the optical wavenumber fco, re- 
spectively. 

Combining Eqs.(|[ [l0|, |ll|), we find the following 
Fourier expansion for the dipole field in the region 
x 3 < z below the dipole: 

E fs (x)= ||^E fs (K)expi(K.X-fc 3a ; 3 ) (13) 

where the 'incident wave vector' is k = (K, — k 3 ) and 
the dipole field's Fourier components equal 

E fs (K) = — !— [fc^d - k(k • d)l expi(-K-R + fc 3 z) 
2e fc 3 

(14) 

Since this vector is perpendicular to k, it may be conve- 
niently expanded into two transverse polarization vec- 
tors e^(K) labelled by the index \i — s, p: 

EfS(K) = ^E e ^( K )MK)-d) (15) 
x expi(— K • R + k 3 z) 

To solve the boundary conditions, the expansion ( p"3| ) 
(and its counterpart ( |l6| ) for the reflected field) is as- 
sumed to be valid down to the surface x 3 = s(X). 
This is actually a hypothesis, the 'Rayleigh hypothe- 
sis', as discussed by Nieto-Vesperinas [[H). Note that 
the present approach also relies on the assumption that 
the dipole is located above the maximum surface height, 
otherwise the absolute value |s(X) — z\ in Eq.(pl|) must 
be retained which complicates the calculation. 

As used by Agarwal [0 we apply the 'extinction 
theorem' |34|] and formulate the boundary conditions as 
integral equations involving the field immediately above 
and below the surface. The integrand is evaluated at 
the surface, and only zeroth and first order terms in 
the profile function s(R) are taken into account. The 
reflected field thus contains zeroth and first order con- 
tributions that are discussed separately in the following. 



3.2.2 Zeroth order: flat surface 

The zeroth order result for the reflected field is obtained 
from the reflection of each Fourier component from a 
flat surface and involves the corresponding Fresnel co- 
efficient r^(K). In the expansion 



E r '°(x) = y|^E r ' (K)expi(K.X + fc 3 z 3 



) (16) 



the Fourier coefficients are thus 



E ,o (K) = i^e^(K)r M (if)(e M (K).dI17) 
3 ^7 
x exp i(— K • R + fc 3 z) 



■5 



where the e^(K) are the unit polarization vectors for 
the specularly reflected waves which are transverse to 
the wave vector (K, k$). From this result we read off the 



reflected Green function for the flat surface and insert 
it into the general formula (^|), giving 



8u9(z) 



1 ST°j(z) 

2 V 



3?ri f d 2 K exp (2ifc 3 z) 



knk 



(18) 



As expected, this result only depends on the dipole's distance z, and not on its lateral coordinate R. The integral 
over the azimuthal angle of the two-dimensional wave vector K may be done analytically, and one finds the familiar 
expressions for a dipole oriented perpendicular (_L) or parallel (||) to the surface [Q: 



SuPjz) 


i 8T\{z) 


3i 
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i STp) 


3i 


Too 


2 Tcxd 
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: exp 



du ■ 



VT 



(2ifc 2\A ~ u 2 ) 



((1 — u 2 )r v — r s ) exp (2ifc z\/ 1 — u 2 ) 



(19) 
(20) 



The integration variable is the reduced wave vector 
u = K/ko. Note that the integration range < u < n 
corresponds to field modes that are plane waves in at 
least one half-space: the lincwidth only depends on 
these modes. For modes with larger wave vectors, 
u > n, the phase factor e 2lk3Z and the reflection co- 
efficients r M become real and the integrands in Eqs.(|l9|, 
p0[ ) become purely imaginary. Hence these modes only 
appear in the lineshift. This property does not hold any 
more when absorption in the dielectric is taken into ac- 



count [[26|, |27|] because the refraction index n and hence 
the reflection coefficients r M are complex for any wave 
vector. 



3.2.3 First order: lateral structure 

The first order contribution to the reflected field is 
Fourier-expanded as in Eq. (|l6|) , and from the calcula- 
tion outlined above, one finds that the Fourier compo- 
nents equal |30|, |35|. |3(| 



E r,1 (K') - i(n 2 - l)L(K') ■ / A^_fc oS (K' - K) E tr '°(K) expi(-K ■ R + k 3 z). (21) 

J (2tt)^ 



We use the notation K for the wave vectors of the zeroth-order field, as in Eq.(|l7|), while K' denotes the wave 
vectors of the scattered field. In Eq.(|2l]), s(K' — K) is the Fourier transform of the surface profile, and E tr '°(K) is 
the Fourier component of the field transmitted by the flat surface given by the Fresnel coefficients t^{K). Finally, 
L(K) is a 3 x 3-matrix given by 

T (TC\ fc " pll , - K * K > - (^3^3, + k^KjSsi) 

L ^> ~ h + k 3 Jv + kQ(n 2 k3 + k3n) W 



where P^- 1 - are projectors parallel and perpendicular to the a;y-plane (the mean surface) and k 3n = \J n 2 k^ — K 2 is 
the vertical component of the transmitted wave vector. It is understood that the third component of the in-plane 
vector K vanishes. 

From the first-order reflected field (|2l| ) we find the following contribution to the Green function: 

- i( ^-l)^/||x: i(l (K',/0J-MK'-K, 

x J2 <t(K)t„(K)e„. : ,(K)expi(K'-X-K-R + t;i3 + t32) (23) 

H=a, p 

where the unit polarization vectors eJ; r (K) describe the bution to the linewidth and the lineshift. 
field transmitted through the flat surface. Upon inser- 
tion in our formula (^|) , one finds the first-order contri- 
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3.3 Transfer function 



We observe in Eq.(23) that each wave vector K of 
the free-space dipole field is diffracted by the Fourier 
component s(K' — K) of the surface profile in such 
a way that the propagation from the dipole down to 
the surface and back again gives rise to a phase factor 
expi[(K' - K) • R + (k' 3 + k 3 )z\. This leads to a lat- 
eral modulation of lineshift and -width at the 'grating 
vector' Q = K' — K. It is expedient to choose this 
wave vector as integration variable in (^3|). One ob- 
tains the following form for the first-order contribution 
to linewidth and -shift: 

Scjj(r) ifl*(r) f d 2 Q 



2 r„ 



/ (^ fc os(Q)^(Q;z)expi(Q-R) 

(24) 



where a dimensionless transfer function 

3Tr(n 2 - 1) f & 2 K expi(/c 3 + k 3 )z 



(25) 



2 J (2tt) 2 k k 3n 
x E E^( K X%( K )We M (K) 

H=s, p k 

has been introduced. In this formula, it is understood 
that K' = K + Q and k' 3 is the corresponding verti- 
cal wave vector component (eq. (p"2|) ) . Eq.(p4|) shows 
that Fj(Q) determines the relative contribution of the 
profile's Fourier components s(Q) to the linewidth and 
-shift. Of particular interest is the width of this "filter" 
as a function of the grating vector Q since it determines 
the lateral resolution of the image. 

We show in Fig.|| contour plots of the imaginary 
part of the integrand in Eq. ( |25| ) for two different grat- 
ing vectors Q. It is apparent that the integrand does 
not have a simple angular dependence in the plane of 
wave vectors K. This implies that in contrast to the 
flat surface case Eq. (|lj) , the angular integral cannot be 
done analytically here, and the transfer function has to 
be computed numerically. In Fig|^ one or two circular 
structures appear to dominate the integrand, depending 
on the magnitude of the grating vector Q. To interpret 
these features, we come back to the Green function ( p3| ) 
that describes the first-order reflected field. For a given 
Fourier component K of the incident dipole field and a 
given grating vector Q = K' — K of the surface profile 
the scattered field amplitude is equal to the profile's 
Fourier amplitude s(Q), multiplied by the product of 
two factors. The first is an electromagnetic scattering 
factor represented by the matrix L(K'), the transmis- 
sion coefficients t^K) and polarization vectors e M (K) 
for the flat surface. The second factor is the exponen- 
tial expi(fc3 -(- k' 3 )z describing the vertical propagation 
of the field from the dipole to the surface and back 
again. The magnitude of the second factor depends on 



the character (propagating or evanescent) of the inci- 
dent and diffracted waves. In particular, the wave vec- 
tors k 3 and k 3 become imaginary for large K, K', and 
the exponential is very small. This limits the relevant 
range of wave vectors that contribute to the integral. 

In the case of the linewidth, the limitation is even 
more strict. It is determined by the imaginary part of 
the integrand (cf. Eq.(|25|)) and is given by two circu- 
lar domains of incident wave vectors K with K < nko 
or K' < nko because it is only in these domains that 
the exponential expi(fc3 + k' 3 )z and the electromagnetic 
scattering factor become complex. These regions are 
clearly visible in the right panel of Fig.||. The grating 
vector Q is here sufficiently large to separate the two 
disks. The disks are merged in the left panel because 
the grating vector is smaller. To summarize, the lateral 
variation of the linewidth above a corrugated surface 
is dominated by two different processes: propagating 
Fourier components of the dipole field are diffracted 
into evanescent waves and interfere with the dipole field 
(the right circular disk of Fig.|], centered at K = 0) and 
conversely, evanescent Fourier components of the dipole 
field are diffracted into propagating waves (the left disk, 
centered at K = — Q). 
3.4 Scanning modes 

Up to now we have determined the fluorescence spec- 
trum of a molecule at a constant height z above the 
structured surface. It is also possible to perform these 
calculations for the more common SNOM scheme of 
'constant-gap' where the separation of the molecule 
from the underlying surface profile is kept at a value 
d. In the notation of the present paper one measures 
the quantity Tj(TL,d + s(R)) where z = d + s(R) is 
the vertical coordinate of the molecule. When calcu- 
lating this linewidth from Eq.(||) one has to take into 
account that our theory only describes surface corruga- 
tions s(R) small compared to the gap d. We thus find 
a linewidth 



r J -(R,d+ s (R)) = r c 



dST°Ad) 

-ff d +» R J + 6T ( R > d ) 

J oz J 

(26) 

The lateral structure is contained in the last two terms, 
the first of which corresponds to the derivative of the 
flat-surface linewidth. It turns out that the transfer 
function for the constant-gap mode can be written in 
the following form 



F j (Q;d)=F j (Q;d)-F j (0;d) 



(27) 



where in the second term the constant-height transfer 
function ( |25| ) is evaluated at zero wave vector. For sim- 
plicity, we focus on the constant-height mode in the rest 
of this paper . 
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Figure 3: Integrand of the transfer function F^Q; z) ( Pq ) as a function of K in units of the optical wave vector ko, 
for fixed grating vector Q||e x . Left panel: grating period comparable to wavelength, Q = O.Sfcoe^. Right panel: 
subwavelength grating period, Q = 5koe x . The contours show the imaginary part, with dark shading indicating 
large negative values. The inset shows a cut through the dotted horizontal line (solid line: negative imaginary part, 
dashed line: real part). 

The molecule's distance from the mean grating surface is kgz = 0.3. Its dipole moment oscillates perpendicular to 
the mean grating surface (along the 23-axis). The substrate index is n = 1.5. We plot the integrand multiplied by 
(2irk ) 2 . 



4 Imaging a grating 

We now examine linewidths and -shifts when the dipole 
is laterally scanned at a constant height above a si- 
nusoidal grating with surface profile s(R) = s(x) = 
hcosQx (grooves parallel to the y-axis). This sim- 
ple geometry reveals the dependence of the radiative 
properties on the four most important length scales (cf. 
fig.|l|): the corrugation height h, the corrugation pe- 
riod a (equal to 2tt/Q for a grating), the molecule's 
distance z from the average surface and the transition 
wavelength A = 27r/fco- Finally, one also has to take into 
account the dipole's orientation. Translational symme- 
try implies that linewidth and -shift are independent of 
y, but they are sinusoidal as a function of the lateral po- 
sition x since they depend linearly on the surface profile 
(cf. Eq.(p^) and Fig.||). All the relevant physics is thus 
encoded in their modulation amplitude |37 . Since this 



amplitude is simply proportional to the grating height 
h, this length scale is already dealt with. 

In FigJ^ the grating has subwavelength corrugation 
amplitude (0.016 A) and period (0.1 A), the substrate is 
a dielectric with refractive index n = 1.5 (glass) and 
negligible absorption, while the dipole is polarized per- 
pendicular to the grating surface (z-polarization). We 
observe that at an average distance from the grating 
of 0.032 A, the linewidth modulation amounts to 20 % 
of the natural linewidth. A much larger modulation is 



observed in the lineshift which amounts up to several 
natural linewidths. We would like to stress this feature 
because frequency shifts have not been considered very 
much in the optics community, perhaps because they 
are more difficult to measure. Our calculations show 
that in near-field optics, the lineshift is much more sen- 
sitive (on an absolute frequency scale) to the surface 
corrugation than the linewidth. This relatively large 
effect is related to the large frequency shift close to a 
flat surface, as discussed in subsection 4.2. We discuss 



the dependence of linewidth and lineshift on the dis- 
tance from the grating and its period in the next two 
sections. 

4.1 The linewidth 

In fig.[|(a) we show the amplitude of the linewidth mod- 
ulation above a grating with subwavelength period, as 
a function of the distance z. One observes that the 
three dipole orientations show different behavior, and 
that the linewidth modulation decreases rapidly with 
increasing distance from the grating. This decrease is 
quite well fitted with an exponential law e~Q z , as shown 
in the inset. Such a law is to be expected since as dis- 
cussed at the end of section |3.3| , the linewidth is domi- 
nated by diffraction processes where the incoming wave 
has a small parallel wave vector (see fig-[| right panel). 
The diffracted wave then has a parallel component with 
wave vector K' w Q whose amplitude decays exponen- 
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Figure 4: The modulations in the linewidth 8T 1 {x 1 z) (a) and lineshift Sui 1 (x,z) (b) of a dipole with vertical po- 
larization at a fixed height above a sinusoidal grating. The grating is made in a glass substrate {n = 1.5) and has 
period a = A/10 and amplitude h = 0.1 X/(2tt) w 0.016 A. Note the difference in scale between linewidth and -shift. 
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Figure 5: The amplitude of the linewidth (a) and lineshift (b) modulations vs. distance from a dielectric sinusoidal 
grating (index n — 1.5) with period a = A/10 and amplitude h — 0.1 A/ (2tt) » 0.016 A. The three linear polarizations 
are shown. The insets give the same data on a logarithmic scale (the absolute value is taken). In (a), the lines are 
fits to simple exponentials e~® z , the amplitude being the only free parameter. In (b), the lines are fits to the model 
function K2{Qz) / z 2 introduced in eq.([33|). For clearness, the x-polarization is omitted from the inset. See text for 
more details. 
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tially oc e~® z if the grating vector Q is much larger than 
fco. This is the phenomenon which allows one to exploit 
the modifications of the lifetime to image nanometric 
structures in the near field. Fig.|5|(a) also shows devi- 
ations from a pure exponential decay of the linewidth, 
we come back to these in eq. (p9[) . 

In Fig.^(a) we plot the amplitude of the linewidth 
modulation as a function of the average distance z for 
a large grating period. One observes that the modu- 
lation amplitude decreases more slowly than in fig.||(a) 
and even shows oscillations for distances of the order of 
the wavelength. In the extreme case of grating period 
much larger than the wavelength this behavior of the 
linewidth may be understood in a simple manner. In 
this case one may write the linewidth modulation in the 
form 

dST° 

a > A : ST(x, z) » 5T°(z-s(x)) « 5T°(z)-s(x)-^- 

(28) 

The linewidth modulation amplitude turns out to be 
proportional to the derivative of the flat surface result. 
In Fig.||(a), the result of this model is indicated by the 
thin lines which coincide quite well with the full calcu- 
lation (symbols) although the grating period is taken 
to be only 5 A. We note, however, that in this simple 
model the two lateral polarizations x and y are always 



degenerate since they both derive from the linewidth 
<5rjj(z) of a dipole polarized parallel to a flat surface. 

Fig.^(a) shows the modulations of the linewidth as 
a function of the grating vector Q for a fixed height z. 
We observe that grating periods larger than A (Q <C fco) 
yield a result similar to that of a flat surface, x- and y- 
polarization being degenerate. As the period decreases 
below the wavelength, this degeneracy is lifted, and we 
observe an overall increase in the linewidth modula- 
tion with some steep features for Q < 2nfco. For still 
smaller grating periods the linewidth modulation de- 
creases again in an approximately exponential manner. 

It is possible to give an asymptotic expression for 
the transfer function ( p5| ) covering the regime of sub- 
wavelength corrugations which is particularly interest- 
ing for applications in optical near-field microscopy. 
This asymptotic expansion is motivated by Fig.|| (right 
panel) where we have seen that the integrand of the 
transfer function is dominated by two circular regions 
of wave vectors. In the limit of Q large compared to fco 
these regions are well separated and have approximately 
circular symmetry. In both cases the angular integra- 
tions may be performed analytically, and one arrives at 
formulae very similar to those for a flat surface (|l^, ^) . 
For the three linear polarizations one obtains: 



Q»fc : Imfj(Q;z)« !^_l e -^ ) I,,,/,.:. 



(29) 



Here, the exponent is ay = 1 for the j = x, z polarizations, and a y = for j = y. If the distance z is much smaller 
than the wavelength, the (complex) dimcnsionlcss functions fj(z) (j — x,y,z) that correct the pure exponential 
decay in eq.d29) are given by ||38| 



KdK , 



I- 2 



fx{z) = 

fv(z) = -(/r-i) 
fz(z) = -4i 



ife3Z 



fc 



K 2 



k 3 + k 3n 2fc (n 2 fc 3 + fc 3n ) 



nko 



K 3 dK fc 3 e ik3Z 
fco n 2 fc 3 + fc 3r , 



nko 



K 3 dK n 2 fc n e ifc3Z 



fco n2 k3 + k 3n 



(30) 
(31) 
(32) 



It can be seen from Fig.^(a) that these asymp- 
totic formulae give an excellent representation for the 
modulation of the linewidth above subwavelength grat- 
ings with periods a < A/5. We may therefore use 
eqs.(p9| |32|) to estimate the lateral resolution SR for 
this type of near-field microscopy: the exponential cut- 
off of high spatial frequencies in Eq.(g2) yields 8R~ z, 
the distance from the surface. Note also that the y- 
polarization gives both a smaller signal and a slightly 
worse resolution compared to the other two polariza- 



tions, which is due to the missing of the factor Q/ko in 
Eq.(29). This is physically plausible because the elec- 
tric field is then continuous across the interface, leading 
to less scattering from the surface corrugation. 

4.2 The lineshift 

Apart from its larger modulation amplitude, the 
lineshift above a sinusoidal grating shows a behavior 
not very different from that of the linewidth. In figj|(b) 
the grating period is subwavelength, and one observes 
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Figure 6: The amplitude of the lincwidth modulation (a) and lineshift modulation (b) vs. distance z from a dielectric 
sinusoidal grating (index n = 1.5) with period a — 5 A and amplitude h — 0.1 A/(27r) ss 0.016 A. Thin lines: simple 
model (|2l) involving the derivative of the flat surface result (see text for details) . 
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Figure 7: Amplitude of the linewidth modulation (a) and lineshift modulation (b) vs. grating vector for a dipole 
at a distance of z = 0.2A/(27r) w 0.032 A from the mean surface. The three linear polarizations are shown, (a) 
The lines give the linewidth modulation obtained from the asymptotic formula (^9j) . (b) The open dots close to the 
thick line give the lineshift for an isotropic dipole, and the thick line itself shows the simple model of Eq.(|33|). The 
electrostatic value of the coefficient C3 is (r oo /8fcQ)(n 2 — l)/{n 2 + 1). 
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again a rapid decrease of the modulation amplitude 
with increasing distance. Vector diffraction is relevant 
and leads to different results for the three polarizations, 
the x-polarization in particular showing a sign change 
(dots). At distances larger than 1/Q, the lineshift shows 
an exponential decrease similar to the linewidth, as 
shown by the solid and dotted lines. These lines are fits 
to the model function K2(Qz)/ z 2 introduced in eq.(|33"|) 
below. 

The case of a grating period larger than the wave- 
length is shown in fig|j](b). The x- and y-polarizations 
show identical lineshifts as above a flat surface. The dis- 
tance dependence is well described by the simple model 
based on eq.(p8|), involving the derivative of the flat- 
surface result. In particular, the lineshift shows a 1/z 4 
power law if the dipole's distance is smaller than about 
A/(27r), as expected from electrostatics. 

Finally, fig.0(b) displays the amplitude of the 
lineshift modulation if the grating period is varied. The 
x, y-polarizations show smooth crossovers from large to 
small periods, and the modulation amplitude globally 
decreases for very small periods. This latter feature can 
be understood from a simple electrostatic calculation, 
as we discuss now. 

We model the substrate as a continuous distribution 
of dipoles (with density Pdip) with which the molecule 
interacts via a (scalar) c^/r 6 law. The total frequency 
shift is obtained by integrating over the half-space 
filled with these dipoles. For a flat substrate one ob- 
tains the familiar power law Soj°(z) = —c^/z 3 with 
C3 = TrpdipCe/6. If the substrate is corrugated, the 
surface region gives an additional contribution to the 
frequency shift. To first order in the surface profile, one 
obtains: 



z<^-: ^(Qjz) 

Z7T 



s{Q)^%K 2 (Qz) (33) 
2 z z 



where <5w 1 (Q; z) is a Fourier transform with respect to 
the lateral coordinates R and Ki(Qz) is the modified 
Bessel function of the second kind. As we show in 
fig.|7](b), this simple model (thick solid line) describes 
quite well the frequency shift for an unpolarized dipole 
(averaged over the three linear polarizations, shown by 
the open circles). For large periods the shift becomes 
independent of Q and tends to —3c3/z 4 , the derivative 
of the flat-surface shift (this follows from the proper- 
ties of the Bessel X-function) while for small periods 
we find an exponential suppression similar to eq . (|29|) : 

Q„ /=r)3/2 -Qz 

Qz » 1 : ^(Q; z) « ~ S (Q) 2% J/2 ( 34 ) 

In this model we may thus explain the suppression of 
high spatial frequencies in the lineshift variations by 



the fact that the frequency shift samples a patch of the 
surface whose radius is of order z, thus washing out 
structures at lateral scales smaller than a < z. As a 
consequence, we expect for lineshift images a lateral 
resolution of the order of z. 

We finally note that similar to the linewidth (p9|), 
the lineshift (34) above a subwavelength grating does 



not show a pure exponential decay with increasing dis- 
tance z (see also inset of fig.|5|(b)). 

5 Imaging an arbitrary substrate 



As pointed out above, our theory is linear in the sur- 
face corrugation and hence able to describe both sinu- 
soidal and arbitrary profiles. An example of a generic 
(two-dimensional) surface is shown in Fig.||. One ob- 
serves a low lateral resolution and quite a weak sig- 
nal at an average height z — 100 nm sa X/2tt (top 
panel) with the linewidth and -shift having compara- 
ble magnitudes. The situation changes dramatically at 
closer distances [z — 20 nm sa 0.2A/27T, middle panel) 
where subwavelength structures are well resolved. Note 
that the linewidth gives a slightly poorer 'image quality' 
than the lineshift. This is due to the fact that the spec- 
tral response of the lineshift behaves more smoothly as 
a function of wave vector than that of the linewidth 
(compare figs.f^(a) and (b)). In other words, some spa- 
tial frequencies are enhanced in the linewidth image, 
leading to a distortion of the observed structures. 

One advantage of optical near-field microscopy over 
other scanning probe techniques is its ability to yield 
information beyond the sample topography, namely 
about its optical contrast. In fact, often samples with 
large topographic features are undesirable in SNOM 
because they lead to the coupling of the optical and 
topographic information |59| ^(J. It is important to 
point out that our theory also applies to substrates with 
purely optical contrast and no topography. For this we 
use the result of Carminati and Greffet that in near- 
field optics variations of the dielectric constant may be 
described by an 'equivalent surface profile' s cq (X) pl[ . 
This quantity corresponds to the vertical integral of the 
optical contrast: 



(X) 



/ 



dx 3 [e(X, x 3 ) - £flat(^3)] (35) 



where £fl a t(£3) is the dielectric function of a flat refer- 
ence substrate. In fig.^J we show an example of such a 
substrate where objects with larger indices are buried 
in a flat substrate. 
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Figure 8: Linewidth and lineshift 'images' of surface topography. The surface profile is shown in the bottom panel 
(thick line, in nm). The mid and top panels show linewidth (thin solid line) and lineshift (dashed line) at two differ- 
ent constant heights above the mean surface, in units of the free-space linewidth r^. Only the laterally modulated 
parts ST 1 (x, z) and 8ui 1 {x, z) are shown. Note the difference in scale for the lineshift at height z — 20 nm. Dielectric 
substrate with index n = 1.5. Transition wavelength A = 628 nm. Vertical polarization. 



6 Limitations of the approxima- 
tion 

The first-order calculation is crucially dependent on the 
linearized boundary conditions at the 'slightly corru- 
gated surface'. More precisely, this means that for all 
relevant Fourier components the following expansion 
must remain sufficiently accurate: 

exp i[(*53 + 4)s(X)] » 1 + i(fc 3 + K)s(X) (36) 

where k 3 and k' 3 are the wave vectors of the incident 
dipole field and the diffracted field, respectively. For 
far-field calculations, k 3 and k' 3 are real and limited to 
the optical wave vector ko. Therefore, we obtain the 
condition |s(X)| < h -C X/2tt, where h characterizes 
the surface corrugation. For calculations in the near 
field one has to include imaginary values of k 3 and k 3 . 
The relevant wave vectors, however, are limited in size 
because the finite distance z leads to an exponential 
damping; this is also apparent in Eq. (p5|) . We hence 



find the condition h -C z. Finally, for a grating with a 
period a well below A all diffraction orders are evanes- 
cent, and the diffracted wave vectors are of the order 
of k 3 ~ Q = 2n/a. In order to perform the lineariza- 
tion d36| ) in this regime, we have to impose the condition 
h <C a. In summary, our method is valid in the regime 

h -C min{a, A, z} (37) 

Note that no restriction is made regarding the relative 
magnitude of the three length scales on the right-hand- 
side. For a sample with optical contrast it is shown in 
Ref.|^l| that the perturbation method is also subject 
to condition (|37|), but now for the equivalent surface 
profile. In particular, this is the case if the index inho- 
mogeneities are confined to a narrow region around the 
interface, below which the sample is homogeneous. 

7 Concluding remarks 

The theory presented here may be generalized to take 
absorption of the substrate into account . In the classical 
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Figure 9: Fluorescence image of a sample with optical contrast. In the dark-shaded region, four objects with a larger 
index n = 2.5 are buried in the flat substrate with n = 1.5. Thin solid line: lincwidth, dashed line: lincshift. The 
other parameters are identical to fig.||. 



picture this is simply done by using a complex index of 
refraction. In the quantum mechanical picture several 
schemes have been proposed fW], to quantize 

the electromagnetic field in the presence of absorbing 
dielectrics, involving different models for the dielectric 
medium. Using the theory of Scheel, Knoll, and Welsch 
(iH), it is easy to check that the fluctuation-dissipation 
theorem still holds. The identification of the field corre- 
lation function in Eq. (Q) with the imaginary part of the 
classical Green function in Eq.(||) hence carries over to 
the absorbing substrate. From the viewpoint of quan- 
tum optics this substantiates the use of the classical 
Lorentz oscillator to compute the fluorescence lifetime 
of real molecules in arbitrary environments. 

We now would like to remark on a few physical ef- 
fects which take place beyond the regime ( p7|) and there- 
fore, are not taken into account by our current treat- 
ment: 

(i) If the grating corrugation is comparable to the 
wavelength h ~ A, one expects many diffraction or- 
ders to be populated. In this regime the calculation of 
the reflected field has to be refined using a full grating 



theory. At large distances from such a 'deep grating' 
the physics should be quite similar to standard far-field 
grating diffraction. At smaller distances the molecule 
samples non-propagating diffraction orders, and evanes- 
cent components of the incidence dipole field could lead 
to qualitative changes of the linewidth. If the grating's 
'depth' exceeds several wavelengths, one may expect 
the formation of a partial photonic band gap, modify- 
ing the near field. Since for a complete band gap only 
evanescent light modes are present, molecular fluores- 
cence would be a very interesting probe to study the 
electromagnetic field in such structures. In a future 
paper we intend to consider a grating with a square 
profile for which an exact diffraction theory is available 

©]47H|H§. 

(ii) If the molecule is put into the selvedge region 
of the grating, i.e. z < h, it is nearly completely sur- 
rounded by the substrate. One then expects that only 
the local environment plays a role, the molecule being 
unable to sense the periodicity of the grating. Numeri- 
cal calculations have been done (l7], [l^, Q which show 
steep variations of the lifetime and a strong polariza- 
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tion dependence. We plan to study the square grating 
model alluded to above to get an analytical insight into 
this situation. 

In conclusion, we have calculated the modification of 
the fluorescence spectrum of a molecule that is scanned 
above a slightly corrugated surface. We have shown 
that the molecule's linewidth and lineshift are influ- 
enced by both surface topography and optical contrast. 
The linewidth acquires variations that amount up to 
20 %, while the lineshift varies over as much as several 
natural linewidths. Furthermore, for an arbitrary sur- 
face profile the lineshift shows a slightly better fidelity 
to the sample structure. We have also presented simple 
models and formulae that allow us to obtain an intuitive 
understanding of the our results for corrugations and 
distances both below and above the molecular transi- 
tion wavelength. Perhaps the most important outcome 
of this paper is that the lateral resolution in the 'flu- 
orescence images' we have obtained is of the order of 
the molecule-surface distance. One would then expect 
to reach a molecular resolution in a novel form of scan- 
ning optical microscopy if the probe molecule could be 
brought nearly in contact with the sample. Recent ex- 
perimental progress in the field of single molecule detec- 
tion, spectroscopy and manipulation give a tantalizing 
hope for the realization of this goal in the near future. 
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